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Executive  Summary  This  report  describes  a  study  motivated  by  the  need  to  replace  legacy 
hydraulic  systems  with  sustainable,  lower-cost,  technology  insertions.  There  have  recently 
been  reports  on  demonstrated  electrical  actuator  systems  that  are  able  to  meet  nominal 
requirements  set  forth  by  the  U.S.  Navy  [1],  and  these  systems  have  or  will  soon  begin  to 
undergo  testing5.  The  goal  of  this  study  was  to  take  a  different  perspective  on  the  retrofit 
problem,  working  instead  on  evaluating  means  for  generating  and/or  verifying  retrofit  or  new 
system  design  requirements.  The  intent  was  to  develop  physics-based  modeling,  simulation, 
and  design  techniques  that  could  support  the  efforts  being  taken  to  develop  and/or  identify 
appropriate  control  surface  actuator  technologies. 

Two  primary  approaches  were  taken.  First,  an  integrated  submarine/actuator  subsystem 
model  was  developed  as  a  platform  for  simulation-based  design  requirements  generation. 
The  utility  of  this  method  was  demonstrated  by  determining  the  power  and  energy  storage 
requirements  for  a  control  surface  actuation  system  by  simulating  known  extreme  subma¬ 
rine  maneuvers,  such  as  emergency  surfacing.  Example  power  requirements  for  an  actuator 
system  for  a  Virginia  class  submarine  were  found  to  be  150  kW  in  four  quadrant  operation 
(with  a  power  limit  of  237  kW)  during  full  sweep  maneuvering.  Additionally,  estimates  of 
local  energy  storage  requirements  in  the  event  of  systematic  power  failure  were  found  to  be 
80  kJ  for  an  emergency  surfacing. 

This  study  also  investigated  synthesis  methods  for  realizing  possible  actuator  system  config¬ 
urations  based  on  performance  specifications.  The  intent  of  this  work  was  to  evaluate  tools 
that  could  support  the  retrofit  process.  The  specific  methods  adopted  are  based  on  classical 
impedance  synthesis  techniques,  first  used  in  design  of  passive  and  active  electrical  circuits. 
In  this  work,  it  is  demonstrated  how  these  concepts  can  be  extended  to  provide  design  candi¬ 
dates  for  electromechanical  systems,  as  might  be  needed  for  control  surface  actuation.  This 
study  demonstrates  how  actuator  system  design  candidates  can  be  formulated  by  consid¬ 
ering  the  case  of  ocean  disturbances  rejection  during  control  of  a  nominal  trajectory.  The 
application  to  disturbance  rejection  shows  how  these  techniques  can  provide  design  alterna¬ 
tives  that  achieve  this  function  either  through  a  purely  passive  design,  which  consumes  no 
additional  power,  or  by  using  a  fully  active  electromechanical  system.  These  methods  can 
provide  alternatives  to  systems  proposed  by  domain  experts,  and  enable  means  for  choosing 
different  configurations  and  initially  sizing  system  components. 


5Reported  February  2011:  http:  / / www.onr.navy.mil/en/Media-Center /Fact-Sheets /Payload-Tube- 

Elect  ric- Act  uat  or .  aspx 
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As  a  current  priority,  the  Navy  is  looking  into  methods  to  reduce  life  cycle  costs  for  future 
vessel  deployment  [1,  2],  One  such  aspect  of  cost  reduction  is  the  replacement  of  existing 
onboard  control  surface  hydraulic  actuator  designs,  which  have  been  in  service  since  World 
War  II,  with  electromechanical  actuators.  Control  surface  actuation  for  shipboard  systems 
has  relied  heavily  on  hydraulic  systems,  taking  advantage  of  legacy  design  knowledge.  How¬ 
ever,  these  hydraulic  systems  prove  costly  to  construct  and  sustain.  Initial  studies  suggest 
whole  life  costs  (WLC)  could  be  reduced  by  up  to  50%  [3]  due  to  conversion  from  hydraulic 
to  electromechanical  systems. 

A  recent  physical  retrofitting  of  a  hydraulic  actuator  on  a  Type  23  Frigate  with  an  elec¬ 
tromechanical  system  was  reported  by  Stafford  and  Osborne  [4],  Two  linear  actuators  were 
used  to  control  the  aft  fixed  fin  stabilizers  on  the  vessel.  This  appears  to  be  one  of  the  first 
instances  of  an  electromechanical  retrofit  on  a  naval  vessel  reported  in  the  open  research 
literature.  The  results  of  the  study  demonstrate  not  only  the  feasibility  of  the  retrofitting 
process  but  also  the  viability  of  electromechanical  systems  in  the  naval  setting. 

Further  justification  for  continued  deployment  of  these  systems  should  show  the  ease  of  main¬ 
tainability  of  electromechanical  systems  over  hydraulics.  Hydraulic  systems  entail  additional 
piping  and  ancillary  power  conversion  systems,  which  add  considerably  to  both  construction 
and  maintenance  overheads.  Hydraulic  fluids  are  drained  and  the  system  is  flushed  in  a 
process  that  takes  upward  of  48  hours  in  addition  to  any  typical  servicing  which  would  be 
standard  for  any  actuation  system  [4] .  Moreover,  modern  electrical  actuators  employ  intelli¬ 
gence  maintenance  monitoring  with  integrated  sensors  and  the  ability  for  prognostics  on  the 

fly- 

Additionally,  retrofitting  with  electromechanical  systems  is  motivated  by  electrification  of 
shipboard  systems,  which  promises  improvements  in  actuation  systems  and  resolution  of 
problems  associated  with  hydraulic  systems.  Not  only  do  hydraulic  systems  add  additional 
energy  conversion  paths  in  the  submarine  system,  additional  crew  training  beyond  electrical 
knowledge  is  required  for  maintenance  and  emergency  repairs.  The  study  by  Stafford  and 
Osborne  [4]  found  that  an  electromechanical  conversion  required  only  two  hours  a  month  for 
planned  maintenance  requirements,  with  other  hourly  reductions  resulting  from  considerable 
time  savings  during  dockside  maintenance  [4],  In  modern  submarines,  actuators  serve  as 
a  synthesis  of  functions:  motion  conversion,  power  conversion,  energy  storage  and  system 
control  as  shown  in  Figure  1.  Due  to  the  hostile  nature  of  the  ocean,  control  surface  actuators 
must  function  in  strenuous  circumstances  and  ensure  functionality  for  surfacing  in  event  of 
power  failure.  Control  surface  actuators  must  also  be  capable  of  meeting  rigorous  torques 
demands  during  diverse  actuation  maneuvers  while  maintaining  high  levels  of  reliability  and 
impact  tolerance.  Other  than  situations  of  overheating  and  extended  over-loading  conditions, 
electrical  actuators  for  control  surface  replacement  attained  98.5  percent  availability  [4], 

Furthermore,  electromechanical  actuators,  which  are  compact  and  lightweight,  are  ideal  can¬ 
didates  to  replace  current  hydraulic  actuators.  Due  to  the  relative  simplicity  and  compact¬ 
ness  of  electromechanical  systems,  the  retrofitted  system  should  yield  a  reduction  in  total 
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Figure  1:  Overall  actuation  system  with  vital  subdivisions. 


actuation  system  weight.  The  proposed  electromechanical  system  tested  by  Stafford  and 
Osborne  resulted  in  a  500  kg  weight  reduction  over  the  existing  hydraulic  system  [4],  The 
conversion  of  electrical  voltage  generated  by  the  submarine’s  power  plant  to  hydraulic  pres¬ 
sure  results  in  energy  loss.  Since  all  energy  conversion  processes  involve  a  loss  in  available 
energy,  removal  of  unnecessary  conversion  processes  results  in  energy  conservation.  One 
potential  cause  of  concern  for  electromechanical  actuation  systems  is  the  continuous  power 
consumption  present  due  to  actuation  position  holding.  In  contrast,  hydraulic  systems  pro¬ 
vide  continuous  position  holding  with  less  power  consumption.  However,  the  experimental 
actuation  fitting  reported  in  [4]  resulted  in  an  electromechanical  actuator  consuming  11.16 
MJ  in  comparison  to  the  standard  hydraulic  actuation  system,  which  consumed  23.76  MJ, 
producing  a  53.03%  reduction  in  energy  consumption  [4], 

Since  hydraulic  actuation  in  submarines  has  an  extensive  history  and  many  system  designs 
and  controls  have  been  determined  using  hydraulics,  a  technique  to  determine  actuator  spec¬ 
ifications  and  component  topology  is  desired.  By  characterization  of  the  existing  hydraulic 
actuators  in  both  time  and  frequency  domains,  electromechanical  actuators  with  identical 
characteristics  can  utilize  existing  control  algorithms  with  indistinguishable  responses.  Addi¬ 
tionally,  investigation  of  the  actuation  systems  must  consider  submarine  system- wide  effects 
on  power  consumption  and  motion  dynamics.  Considerable  concentration  of  mass  towards 
the  aft  section  of  the  ship  results  in  stability  issues  for  submarine  motion. 

The  Stafford  report  indicates  strong  evidence  for  the  favorability  of  electromechanical  system 
retrofitting.  The  report  presents  experimental  confirmation  of  the  possibility  for  a  retrofit 
of  existing  hydraulic  systems  with  positive  benefits.  However,  many  issues  and  concerns  re¬ 
main  unanswered  for  a  finalized  retrofit  solution  such  as  system  component  topology,  energy 
storage  requirements/placement  and  power  consumption  during  intense  maneuvers.  This 
report  addresses  each  of  these  problems  in  turn  and  explores  potential  solutions  for  a  final 
retrofit  design  in  accordance  to  the  specified  parameters  in  Table  1.  Since  the  retrofitting  of 
the  control  surface  actuators  suggests  system  insertion,  there  is  a  need  for  methods  to  aid 
formulation  of  requirements  for  actuation  subsystems  that  correspond  to  submarine  opera¬ 
tional  requirements.  Modern  modeling  and  simulation  tools  enable  a  total  systems  approach 
that  allows  the  possibility  to  propagate  overall  system  response  requirements  into  specified 
subsystems.  However,  early  or  concept  design  applications  requires  that  physics-based  mod¬ 
els  be  developed  with  minimal  design  information.  Detailed  selection  and  sizing  require  that 
requirements  be  formulated,  and  it  can  be  helpful  to  have  systematic  ways  to  model  and 
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Table  1:  Specified  Performance  Parameters  [1] 


Parameter 

Specification 

Operating  Torque  Capacity 

700,  000  ft-lbs  (950  kN-rn) 

Peak  Torque  Capacity 

1,000,000  ft-lbs  (1,360  kN-rn) 

Operational  Angular  Displacement 

±35  degrees 

Rate  of  Angular  Displacement 

5  ±  1  deg/s 

Position  Accuracy 

±0.25  degrees 

Energy  Storage 

1.4  MJ 

simulate  systems  for  gaining  insight  during  this  process.  The  following  section  describes  the 
physics-based  model  for  the  submarine  model  system,  which  forms  a  basis  for  generating  key 
performance  requirements  from  first  principles. 


2  Submarine  Model  Development 


Physics-based  models  can  be  helpful  in  actuator  design  assessment  since  they  enable  critical 
power  levels  to  be  estimated  through  simulation  studies.  This  can  be  especially  important 
in  actuation  and  energy  storage  sizing,  since  designers  may  need  to  consider  power  demands 
imposed  by  a  broad  range  of  external  effects  [5].  Designers  can  virtually  experiment  by 
imposing  changes  in  the  system  components  and/or  environmental  conditions  as  they  see 
fit.  In  actuator  retrofit  studies,  the  dynamics  of  an  existing  system  (e.g.,  a  legacy  hydraulic 
actuator  system)  may  rely  on  a  power  distribution  that  is  radically  different  from  the  novel 
architectures  being  considered,  which  could  mask  new  design  requirements  [4],  For  instance, 
the  source  impedance  of  a  dc  bus  system  remains  intrinsically  different  to  a  hydraulic  supply 
line  [6].  Physics-based  models  can  help  reduce  or  eliminate  any  inherent  bias  or  uncertainty 
associated  with  this  type  of  retrofit  problem.  Furthermore,  as  there  may  be  little  to  no 
information  available  about  how  an  existing  system  was  specified,  a  physics-based  model 
can  provide  a  more  objective  basis  for  making  future  design  decisions,  including  a  way  to 
derive  possible  duty  cycle  descriptions  for  onboard  actuator  systems.  The  latter  would  offer 
an  alternative  or  basis  for  comparison  to  costly  full-scale  sea  trials. 

This  section  summarizes  the  submarine  and  control  surface  actuator  models  developed  for 
the  purpose  of  this  retrofit  study.  The  simulations  are  based  on  high  fidelity  physics-based 
models  meant  to  provide  insight  into  the  interaction  between  the  submarine  and  actuator 
dynamics.  The  simulations  can  be  used  to  guide  electromechanical  actuator  retrofit  require¬ 
ments.  However,  these  results  are  only  a  first  step  toward  a  retrofit  design,  the  next  step 
being  methods  that  could  support  design  of  a  specific  actuator  system,  such  as  described  in 
Section  4.2. 
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Earth-Fixed 


p  -  roll 


Figure  2:  Coordinate  frame  definitions  of  submarine  vehicular  dynamics. 


2.1  Six  Degree  of  Freedom  Submarine  Model 

A  detailed  model  of  a  full  six  degree  of  freedom  (DOF)  submarine  system  was  developed 
and  a  simulation  implemented  using  the  MATLAB/Simulink  environment6.  For  the  char¬ 
acterization  of  energy  requirements  of  actuators,  all  six  degrees  of  freedom  were  included 
to  account  for  all  possible  internal/external  power  flow  requirements.  The  body-fixed  and 
global  reference  frames  used  in  the  derivation  are  illustrated  in  Figure  2.  The  high  fidelity  is 
necessary  as  model  simplifications  will  result  in  gross  underestimates  of  power  requirements. 
As  a  result,  domain-specific  knowledge  of  the  system  is  mandatory  for  an  accurate  design 
assessment.  A  detailed  discussion  of  the  derivation  is  provided  in  Appendix  A.  The  force 
and  moment  models  described  in  Appendix  A.l  form  the  total  forces  and  moments  exerted 
on  the  submarine  body.  A  mathematical  model  in  the  form  of  nonlinear  ordinary  differential 
equations  of  motion,  with  respect  to  a  body-fixed  coordinate  system,  for  the  submarine  takes 
the  form, 

Mv  +  C{v)v  +  D(u)v  +  g{r])  +  g(y)ui  =  t  (1) 

where  Mv  is  the  total  inertial  matrix,  C(u)  is  the  Coriolis/centripetal  matrix,  D{v)  is  the 
total  damping  matrix,  g{ifj  are  the  buoyancy  and  weight  components,  g{v)  is  the  control 
input  matrix,  and  r  is  external  disturbances.  These  equations  summarize  the  external  forces 
imposed  on  the  actuation  units.  The  following  subsections  describe  how  external  forces 
transfer  to  the  internal  actuation  systems. 


6The  Mathworks,  Inc.,  3  Apple  Hill  Drive,  Natick,  MA. 
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Figure  3:  Hydrodynamic  moment  unit  definition. 
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Figure  4:  Hydrodynamic  lift  force  and  moments  imposed  on  control  surface. 


2.2  Control  Surface  Hydrodynamics 

Models  for  representing  the  control  surface  hydrodynamic  interactions  and  actuation  system 
were  also  formulated.  The  importance  of  the  model  equations  is  twofold.  First,  the  control 
forces  imposed  on  the  submarine  motion  can  be  determined,  see  Figure  3.  Second  and  more 
importantly,  the  reaction  forces  imposed  on  the  actuators  due  to  submarine  motion  and 
ocean  disturbances  can  be  found.  These  force  transmissions  are  key  for  the  identification 
of  actuator  power  requirements  [5].  Much  akin  to  an  airfoil,  the  control  surface  experiences 
both  hydrodynamic  lift  and  drag,  see  Figure  4.  Since  we  operate  under  the  assumption  that 
the  control  surface  remains  submerged  for  the  duration  of  the  simulation,  modeling  of  the 
hydrodynamics  was  implemented  with  nonlinear  lift  and  drag  coefficients.  A  more  detailed 
discussion  of  the  control  surface  hydrodynamics  can  be  found  in  Appendix  A. 2.  In  addition 
to  the  hydrodynamic  reaction  forces  discussed  above,  ocean  wave-induced  disturbances  can 
inject /remove  additional  energy  from  the  control  surface  system.  For  the  purposes  of  this 
study,  ocean  disturbances  were  modeled  using  a  Pierson-Moskowitz  spectrum.  This  model 
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Figure  5:  Frequency  response  due  to  Pierson-Moskowitz  ocean  spectrum. 
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Figure  6:  Simplified  layout  of  control  surface  system. 


was  transformed  via  a  Beaufort  number  relationship  [7]  to  derive  the  spectrum  shown  in 
Figure  5.  This  spectrum  was  used  throughout  this  study  to  model  ocean  disturbances.  The 
methodology  used  can  be  applied  to  explore  the  influence  of  other  sea  states,  besides  that 
defined  by  a  Pierson-Moskowitz  spectrum. 


2.3  Control  Surface  Actuation-Ship  Interconnection 

The  interconnection  or  coupling  between  the  control  surface  forces  and  the  main  ship  motion 
requires  that  a  specific  architecture  be  chosen.  The  mechanical  arrangements  for  submarine 
control  surface  controls  can  take  several  different  forms  (e.g.,  as  described  in  [1]).  For  the 
purposes  of  this  study,  the  somewhat  idealized  layout  illustrated  in  Figure  6  was  used  and 
enabled  us  to  experiment  with  key  dynamic  effects.  Specifically,  in  addition  to  the  external 
loading  from  hydrodynamic  forces,  the  layout  chosen  incorporates  nominal  elastic,  inertia, 
and  damping  effects.  No  specific  actuator  technology  appears  in  the  simplified  layout  of 
Figure  6  since  a  goal  of  this  study  was  to  formulate  requirements  for  actuator  systems 
that  could  be  inserted  in  this  space  to  achieve  overall  submarine  motion  requirements.  The 
actuator  system  layout  was  modeled  using  bond  graph  methods,  which  are  useful  in  modeling 
a  very  wide  range  of  engineering  system  types  and  complexity  [8].  The  bond  graph  model  of 
the  actuator  layout  was  useful  in  two  ways.  First,  the  model  was  readily  integrated  with  the 
overall  submarine  model  so  that  full  system  simulations  could  be  conducted  (as  described 
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Control  surface  transient  deflection  Actuator  Power  Consumption  (kW) 


Control  surface  moment  during  emergency  surfacing  Actuator  Energy  Consumption  (kJ) 


Figure  7:  Dynamic  demands  during  full  sweep  simulation. 


later).  Second,  the  representation  of  the  actuator  system  layout  in  bond  graph  form  also 
facilitated  application  of  impedance  methods  as  presented  in  Appendix  C.2.  The  following 
sections  describe  how  these  model-based  studies  were  conducted  and  used. 


3  Formulating  Requirements  for  Actuation  Subsystems 


Two  techniques  useful  in  retrofit  design  were  explored  in  this  study.  First,  it  was  shown 
how  simulations  of  the  full  submarine  can  be  used  to  derive  required  control  surface  actuator 
forces.  For  specified  operational  scenarios,  different  types  of  requirements  can  be  derived. 
The  example  of  finding  energy  requirements  for  extreme  maneuvers  is  described.  Secondly, 
it  is  shown  how  these  models  can  be  used  to  derive  requirements  in  the  form  of  impedance 
functions  to  aid  design  synthesis. 


3.1  Energy  Requirements  during  Extreme  Maneuvering 

For  a  typical  submarine  mission,  the  actuation  systems  encounter  the  greatest  hydrodynamic 
torque  demands  during  sharp  yaw  or  pitch  banks  [9].  To  demonstrate  prediction  of  the 
dynamics  of  a  submarine  performing  extreme  pitch  bank  maneuvers,  a  full  sweep  and  hold 
of  the  stern  plane  actuation  system  was  simulated.  This  maneuver  is  analogous  to  a  rapid 
rising  and  falling  of  the  submarine  during  operation  without  a  change  in  ballast  pressure. 
Additionally,  the  ONR  utilizes  this  basic  test  for  the  quantification  of  actuation  systems  [1] . 
Figure  7  displays  the  results  of  the  full  sweep  maneuver.  During  the  holding  maneuver  in 
Figure  8,  the  torque  demand  continues.  From  the  persistent  demand  for  torque  actuation,  the 
final  actuation  design  must  be  capable  of  four  quadrant  operation,  meaning  both  actuation 
and  holding.  The  other  simulation  of  interest  is  the  emergency  surfacing  maneuver  [1], 
Unlike  the  continuous  sweep  simulation  above,  the  control  surfaces  during  the  maneuver 
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Submarine  States  during  Full  Control  Surface  Sweep 
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Figure  8:  Selected  submarine  states  for  full  sweep  maneuver. 


perform  precision  tracking  to  orientate  the  submarine  to  the  proper  global  angle  (6)  so  that 
the  vessel  exits  the  water  at  the  proper  angle  for  structural  integrity.  Additionally,  with  the 
large  change  of  buoyancy  after  blowing  the  ballast  chamber,  the  stern  plane  actuators  must 
ensure  that  the  ascent  is  controlled  rather  than  allowing  the  submarine  system  to  destabilize. 
In  the  worst  case,  an  uncontrolled  ascension  could  exit  the  water  perpendicularly  causing 
bodily  harm  to  the  crew  and  severe  damage  to  the  structural  systems.  Selected  states  of 
the  submarine  model  during  the  surfacing  maneuver  are  shown  in  Figure  9.  The  submarine 
model  states  correspond  closely  with  previous  submarine  modeling  studies,  such  as  [5].  Power 
flow  analysis  during  the  emergency  surfacing  maneuver  helps  quantify  power  and  energy 
consumption  requirements.  Results  from  this  type  of  simulation  study  are  summarized  in 
Figure  10. 


3.2  Disturbance  Rejection  Requirements 

Simulation  of  extreme  maneuvers  provides  information  on  the  requirements  for  both  the 
power  demands  of  the  actuator  and  the  local  energy  requirements  for  various  maneuvers. 
These  results  can  be  used  directly  by  domain  experts  to  build  specifications  for  new  actua¬ 
tor  systems.  For  a  synthesis  of  design  candidates  from  simulation  data,  other  techniques  can 
be  applied  [10].  To  illustrate,  another  simulation  study  was  constructed  examining  rejection 
of  ocean  disturbances  by  using  the  control  surface  actuators.  For  this  case,  the  actuator 
package  dynamics  were  studied  only  locally  (i.e.,  the  full  submarine  model  was  not  needed) 
with  imposed  hydrodynamic  loads.  The  results  derived  are  the  nominal  system  response 
characteristics  shown  in  Figure  11.  The  two-port  detailed  derivation  of  this  frequency  do¬ 
main  model,  including  hydrodynamic  Hoerner  linearization,  is  provided  in  Appendix  A. 2. 
Superimposing  the  Pierson-Moskowitz  ocean  disturbance  spectrum  (see  Figure  12)  illustrates 
the  hypothetical  energy  cost  of  rejecting  disturbances  with  an  actuator  system  that  is  purely 
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Submarine  States  during  Emergency  Surfacing 
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Figure  9:  Selected  submarine  states  for  rapid  ascension  maneuver. 
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Figure  10:  Control  surface  deflection  and  dynamic  demands  during  emergency  surfacing. 
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Figure  11:  Command-hold  frequency  response  of  control  surface  system. 


Actuator  Sensitivity  to  Ocean  Disturbance  -  (  a>cs/Td) 


Figure  12:  Pierson-Moskowitz  normalized  to  the  Control  surface  System  Response. 


active;  i.e.,  the  disturbance  forces  are  counter-acted  entirely  by  the  powered  actuator.  As 
shown,  the  magnitude  of  the  response  remains  high  in  this  frequency  band,  so  energy  is 
sapped  from  the  system  and  power  required  for  disturbance  rejection  increases.  It  will  be  il¬ 
lustrated  in  section  4.2  that  a  two-port  impedance  synthesis  method  allows  us  to  recommend 
an  alternative  actuator  system  design  that  can  passively  reject  these  disturbances.  As  such, 
these  methods  provide  a  designer  with  alternative  options  for  addressing  the  various  require¬ 
ments,  such  as  disturbance  rejection,  that  may  need  to  be  considered  during  the  actuator 
system  retrofit  process. 


4  Retrofit  Design  via  Specified  Requirements 

The  following  is  a  demonstration  of  ways  for  generating  design  candidates  given  system 
requirements  derived  from  the  simulation  models  previously  discussed. 


Retrofit  Design  of  Submarine  Actuation  Systems 


Page  13 


Ragone  chart  -  Energy  vs.  Power  delivery  profile 
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Figure  13:  Comparison  of  energy  density  to  power  density  of  energy  storage  devices  [12]. 


4.1  Ragone  chart  design  of  energy  storage 

With  the  estimated  energy  requirements  for  extreme  maneuvers  from  Section  3.1,  energy 
storage  devices  can  be  sized  accordingly  [4,  11],  One  commonly  used  methodology,  albeit  ad 
hoc,  is  technology  identification  through  a  Ragone  chart.  For  storage  of  energy  needed  to 
power  an  electromechanical  actuation  system,  many  suitable  technologies  exist  [12].  Figure 
13  shows  many  of  the  potential  energy  storage  technologies  that  could  provide  emergency 
power  and  regulate  the  voltages  provided  to  the  system  electronics.  The  energy  consumption 
predicted  via  simulation  for  the  single  actuator  can  support  systematic  selection  and  sizing 
of  energy  storage  requirements.  A  summary  of  the  key  simulation  results  is  given  in  Table 
2.  In  the  event  of  power  failure,  additional  energy  may  also  be  needed  in  case  actuations 
are  incomplete.  An  askew  control  surface,  for  example,  could  potentially  drive  the  entire 
submarine  system  into  a  state  of  instability.  For  a  single  actuator,  the  energy  storage  should 
be  five  times  the  encountered  peak  for  an  extension  to  the  amount  of  time  a  maneuver  could 
persist  [13].  The  resulting  predicted  energy  storage  required  is  2465  kJ.  Given  the  simulation 


Table  2:  Fundamental  Simulation  Results  for  Single  Actuator 


Full  Sweep 

Peak  Power 

130  kW 

Total  Energy 

493  kJ 

Emergency  Surfacing 

Peak  Power 

10  kW 

Total  Energy 

80  kJ 

results,  ad  hoc  energy  storage  and  actuator  system  sizing  is  possible  [12].  Identification  of 
viable  storage  technologies  can  be  accomplished  using  the  Ragone  chart  in  Figure  13.  The 
energy  storage  mechanism  must  contain  a  large  energy  density  in  comparison  to  the  relative 
power  density  of  the  system.  These  are  preliminary  results,  and  further  investigation  into 
the  preferred  energy  storage  for  submarine  applications  should  be  conducted. 
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Figure  14:  Flowchart  of  synthesis  procedure. 


The  placement  of  the  energy  storage  device  in  the  system  topology  also  remains  heuristical 
[4,  12].  A  potential  design  could  have  the  energy  storage  in  a  centralized  location  much  akin  to 
the  current  hydraulic  accumulator  designs  [9].  However,  centralized  energy  storage  requires 
power  transport  and  distribution  that  add  losses  to  the  system  [14].  Another  possibility  is  to 
have  localized  energy  storage  for  each  actuation  system.  This  setup  provides  quick  response 
times  but  once  the  power  in  one  area  is  exhausted,  energy  from  other  underutilized  actuators 
proves  difficult  to  relocate  quickly  to  drained  systems  in  the  event  of  power  failure  [15].  Also 
note  that  the  voltage  type  of  the  ship  power  grid  is  also  a  crucial  factor  [1], 

While  the  ad  hoc  design  methodology  produces  feasible  designs,  the  topology  and  energy 
storage  mechanisms  are  not  necessarily  ideal  for  the  submarine  system.  Completion  of  the 
design  process  requires  expert  and/or  a  posteriori  knowledge  of  energy  systems  as  well  as 
the  submarine  specific  application.  The  system  simulations  provide  insight  into  actuation 
energy  requirements  for  physical  submarine  maneuvers  which  aid  in  energy  storage  sizing. 


4.2  Retrofit  Actuator  Design  via  Impedance  Methods 

In  the  following,  the  synthesis  approach  is  shown  to  produce  not  only  the  required  energy  stor¬ 
age  for  a  given  actuation  task,  but  also  potential  designs  for  the  energy  storage  mechanism. 
Rather  than  heuristical  placement  of  components  in  a  given  system,  a  direct  identification  of 
both  component  location  and  size  is  found  [16,  17].  Two-port  impedance  synthesis  requires 
less  direct  involvement  from  system  experts  and  allows  for  novel  design  candidate  generation 
[18].  The  steps  in  this  procedure  are  summarized  in  Figure  14.  For  retrofit  design  generation 
in  the  case  of  the  submarine  actuation  system,  we  reformulate  the  notional  multi-domain 
(bond  graph)  model  to  incorporate  a  two-port  black  box  (see  Figure  15).  With  the  goal  of 
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Figure  15:  Two-port  modification  of  notional  system. 


Desired  system  ( w  cs/Td )  frequency  response 


Figure  16:  Adjusted  desired  response  of  system. 


rejecting  disturbances  having  frequency  content  described  by  a  Pierson-Moskowitz  spectrum, 
a  desired  hold-command  response  can  be  adjusted.  Attenuation  of  disturbance  frequencies 
in  the  band  of  the  Pierson-Moskowitz  spectrum  becomes  a  new  desired  system  response,  as 
seen  in  Figure  16.  Executing  the  synthesis  methodology  with  the  given  system  and  desired 
response  produces  a  two-port  impedance  function,  Z(s)  [10].  Converting  Z(s)  to  a  real  sys¬ 
tem  can  prove  to  be  problematic,  as  discussed  in  Appendix  C.2.  Typically  the  form  of  the 
Z(s)  function  makes  it  difficult  to  formulate  a  physical  realization,  especially  in  the  face  of 
uncertainty  [19].  Additionally,  the  resulting  impedance  function,  Z(s),  contains  both  active 
components  and  passive  components.  The  active  components  of  Z(s)  can  be  implemented 
through  the  actuator  as  feedback  as  illustrated  in  Figure  17.  However,  order  reduction  (bal¬ 
anced  realization)  on  the  system  can  lead  to  a  good  approximation  of  the  system  as  shown 
in  Figure  18.  This  approximate  function  form  of  Z(s)  can  lead  to  a  more  realistic,  simplified 
system  realization.  The  order  reduction  in  this  case  removes  the  need  for  high  order  active 
components  [19].  As  a  result,  the  actuation  system  can  reject  ocean  disturbances  passively. 
The  basic  structure  of  this  design  is  shown  in  Figure  19.  Interestingly,  the  result  of  using  an 
impedance  synthesis  methodology  is  a  design  which  contains  a  local  energy  storage  device. 
This  energy  storage  mechanism  recharges  using  power  from  ocean  disturbances.  The  system 
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Figure  17:  Block  diagram  representation  of  Z(s). 


Figure  18:  Model  order  reduction  for  design  realization. 
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Figure  19:  Block  diagram  representation  of  purely  passive  ocean  rejection  system. 
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Figure  20:  (a)  Linear  mechanical  realization  of  reduced  impedance  function,  (b)  One  poten¬ 
tial  rotational  realization  of  the  reduced  impedance. 


Ky  bi 2  Kjmi2 


rejects  external  energy  by  storing  it  locally  in  a  notional  capacitive  element.  As  shown  by 
Figure  20,  this  energy  storage  mechanism  can  be  realized  in  alternative  energy  domains;  i.e., 
using  different  technological  insertions.  In  the  case  of  this  design  process,  the  push  towards 
an  integrated  electrical  power  system  may  suggest  that  Figure  21  would  produce  a  preferred 
storage  mechanism. 


5  Conclusions 


Both  the  simulation-based  and  impedance-based  synthesis  techniques  developed  in  this  work 
are  viable  approaches  for  aiding  early-stage  design  and/or  retrofit  of  actuator  subsystems. 
The  full  system  model  simulations  can  be  used  to  study  various  operational  maneuvers  so 
that  actuation  system  requirements  can  be  formulated  directly.  Synthesis  techniques  can 
employ  relatively  basic  system  models  and  can  be  useful,  especially  in  light  of  a  limited 
knowledge  base,  in  providing  alternative  design  configurations.  Both  approaches  and  results 
are  reported  in  more  detail  in  [10,  19,  15,  14]. 


In  the  applications  examined,  the  synthesis  approaches  lead  to  compensating  effects  for 
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actuators.  Nevertheless,  it  was  demonstrated  that  the  resulting  designs  could  lead  to  a 
notable  reduction  in  energy  consumption  for  systems  with  identical  power  density.  This 
suggests  that  the  synthesis  methods  can  help  guide  how  passive  and/or  active  topologies  for 
energy  storage  might  be  combined  when  designing  and/or  retrofitting  actuation  subsystems. 
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A  Modeling  and  Simulation  of  Submarine  System 


In  this  appendix,  the  time-invariant  nonlinear  submarine  vehicular  kinematics  and  dynamics 
are  introduced.  For  the  sake  of  the  reader,  the  equations  are  partially  derived  to  reveal  the 
framework  of  the  submarine  motion  and  forces.  Additionally,  a  complete  presentation  of  all 
of  the  equations  with  definitions  of  every  term  can  be  found  in  Gertler  [38]  and  Feldman  [33] 
and  a  comprehensive  derivation  of  all  of  the  equations  can  be  found  in  Watt  [5,  83].  The 
notation  of  derivation  follows  that  of  Fossen  [34]  with  some  modifications. 


A.l  Mathematical  Modeling  of  Submarine  Vehicular  Dynamics 


By  establishing  simple  coordinate  frame  relationships,  the  overall  complexity  of  a  six  degree 
of  freedom  (DOF)  model  can  be  reduced  to  the  least  complex  nonlinear  form.  The  position 
and  orientation  of  the  submarine  are  relative  to  the  earth-fixed  frame  and  the  linear  and 
angular  velocities  are  relative  to  the  body-fixed  frame.  Representing  these  expressions  in 
terms  of  vectors,  the  system  state  variables  are  77  and  v  and  the  forces  and  moments  are  r 
as  illustrated  by  equation  (2). 

rj  —  [  x  y  z  (f  9  i/j  ]T 

v  —  [uvwpqr]T  (2) 

t=[  X  Y  Z  K  M  N  ]T 

The  velocity  vector  of  the  earth-fixed  frame  and  the  velocities  of  the  submarine  fixed  reference 
frame  are  related  by  equation  (3). 

V  =  J(i  i)v  (3) 

where  J (77)  is  the  rotation  matrix  expressed  as  equation  (4). 


J(v ) 


Jl(V 2)  03x3 

03x3  Afe)] 


(4) 


The  subordinate  matrices  J\  (r/2)  and  A(h 2)  are  given  by  equation  (5)  and  equation  (6), 
respectively. 


AM 


cos  i[>  cos  9  Ji  2(7/2)  A3M 

sin  ip  cos  9  J22M  A 3(772) 

—  sin  9  cos  9  sin  cos  9  cos  tp 


(5) 


AM 


1  sin  < p  tan  9  cos  cp  tan  9 

0  cos  <p  —  sin  ip 

0  sin  79/cos  9  cos  92/cos  9 


(6) 
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where  the  remaining  components  of  matrix  (5)  are  given  by  equation  (7). 

J 12  (772)  =  —  sin  ?/?  cos  p  +  cos  7/1  sin  9  sin  f 
J13 (772)  =  sin  7/1  sin  f  +  cos  7/1  cos  p  sin  9 
J22  (7/2)  =  cos  7 1)  cos  f  +  sin  p  sin  9  sin  7/1 
J23  (7/2)  =  —  cos  7/1  sin  p  +  sin  9  sin  7/?  cos  p 


(7) 


A.  1.1  Nonlinear  Dynamics 


Combining  terms  from  the  expansion  of  the  Lagrange  equations,  the  rigid-body  dynamics  of 
the  submarine  are  denoted  by  equation  (8). 

Mbu  +  Cb{v)v  =  t  (8) 


where  the  vector  v  represents  how  the  equation  is  defined  with  respect  to  the  body.  The 
rigid-body  inertia  matrix,  Mb,  is  defined  in  equation  (9).  Cb  represents  the  centripetal  and 
Coriolis  portions  of  the  inertial  matrix  displayed  in  equation  (10). 


m 

0 

0 

0 

mzG 

—myG 

0 

m 

0 

— mzG 

0 

mxG 

0 

0 

m 

myG 

—mxG 
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0 

— mzG 

myG 

h 

—  Ixy 

4Z 

mzG 

0 

—mxG 

—lyx 

4 

—lyz 

~myG 

mxG 

0 

I zx 

—  Izy 

4 

where  m  is  the  mass  of  the  submarine,  /  is  moment  of  inertia  (or  correspondingly,  the  product 
of  inertia)  with  respect  to  a  rotation  axis,  xG  is  the  x-coordinate  of  the  center  of  gravity, 
7 ig  is  the  y-coordinate  of  the  center  of  gravity,  and  zG  is  the  z-coordinate  of  the  center  of 
gravity. 


cb  M  = 


cb2(v) 


cb(u)  = 

m(yGq  +  zGr) 

-m(yGp  +  w) 

— m(zGp  —  v ) 

0 

fyzQ  7zr 
_  -IyzT  +  Iyq 


03x3  Cbi(u) 

cb2(u) 

(10) 

— m(xGq  —  w) 

—m(xGr  +  v ) 

m{zGr  +  xGp ) 

—rn  (yGr  —  u ) 

(11) 

- m(zGq  +  77) 

m(xGp  +  yGq )  _ 

-Iyzq  +  Izr 

1 

1 

is? 

0 

Ixz T  “1“  lx V 

(12) 

fr ::  7  Ixp 

0 

Due  to  the  components  of  vector  v  in  the  Coriolis/centripetal  term  which  couple  with  other 
terms  in  the  v  vector,  the  dynamics  of  the  submarine  are  nonlinear.  Also  note  that  the  Cb 
matrix  is  in  skew-symmetric  form. 
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Since  the  buoyancy  and  the  weight  forces  act  due  to  the  orientation  of  the  submarine  to  the 
earth-fixed  frame,  a  transformation  between  the  earth-fixed  frame  and  the  body  frame  must 
be  preformed.  Using  the  transformation  ,J\  (r/2),  the  restoring  moments/forces  coefficients 
are  given  by  equation  (13). 

(W-B)  sin  6 
—(W  —  B )  cos  6  sin  p> 

—  (W  —  B)  cos  6  cos  <p 

—  ( UgW  —  UbB)  cos  6  cos  ip  +  (zqW  —  zbB )  cos  6  sin  <p 
[zqW  —  ZbB )  sin  0  +  (xqW  —  xbB)  cos  6  cos  p> 

—  {xqW  —  XbB)  cos  6  sin  >p  +  {ycW  —  ysB)  sin  0 

where  (xq,  ye,  zq)  and  (xb,  yB ,  zb)  are  the  center  of  gravity  and  center  of  buoyancy, 
respectively. 


The  potential  damping  coefficients  clue  to  the  skin  friction  are  found  as  in  Watt  [83]  as  the 
integral  of  the  viscous  stress  tangent  to  the  submarine  body.  As  in  the  derivation  of  Watt  [83] , 
the  terms  higher  than  second  order  are  negligible  and  can  thus  be  dropped.  Additionally,  clue 
to  the  planes  of  symmetry  of  the  submarine  body,  the  damping  matrices  may  be  simplified 


to, 


D{y) 


Dn(v)  D12(u) 
03x3  D22(v) 


(14) 


Du(v) 


Xu  T  AuH  \u\  0  0 

0  Yv  +  Yv |q  |u|  0 

0  0  -|-  Yu!ju!j  |  w 


(15) 


0  Xq\q |  |g|  0 

D12(u)  =  Yp\p\  \p\  0  Yr |r|  \r\  (16) 

0  Yq |,|  |g|  0 


I\.p  T  hp\p\  \p\  0  Rr\r\  |r| 

D22{y)  =  0  Mq  +  Mq\q\  |g|  0  (17) 

Yp\p\  \p\  0  Nr  +  iVr|r|  \r\  _ 


where  the  first  terms  in  the  matrix  are  the  coefficients  of  the  skin  friction  factors  and  the 
second  terms  are  the  coefficients  of  the  vortex  shedding  drag  factors. 


The  force  components  of  the  control  inputs  can  be  expressed  as  a  product  of  the  input  matrix 
and  the  input  vector, 

Ti  =  g{v)ui  (18) 

where  g{y)  is  the  control  parameter  matrix  which  contains  the  parametric  relationships  of 
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the  input  vector  rq.  The  control  parameter  matrix  g{y)  is, 


9(v) 


0  0  0  1 

Yuu5ru2  0  0  0 

0  Zuus.u  A  u  uSif  ^  0 

0  0  0  0 

0  MuuSsu 2  MuuSbu 2  0 

Nuusru2  0  0  0 


(19) 


where  the  control  parameters  are  given  by  the  size  and  shape 
input  vector,  rq,  is  given  as, 


u,  = 


5r 

5S 

5b 

v 

|_  -S'- prop 


of  the  control  surfaces. 


The 


(20) 


where  5r  is  the  deflection  of  the  rudder  in  radians,  5S  is  the  deflection  of  the  stern  plane  in 
radians,  5b  is  the  deflection  of  the  bow  plane  in  radians  and  Xprop  is  the  simplified  thrust 
model  which  only  applies  force  in  the  surge  direction. 


A. 1.3  Full  Equations  of  Motion 

By  compiling  the  forces  and  moments  derived  in  the  preceding  analysis,  the  total  forces 
and  moments  on  the  submarine  body  can  be  found.  The  total  inertial  matrix  and  Cori¬ 
olis/centripetal  matrix  can  be  found  by  summing  both  the  body  dynamics  and  the  added 
mass  matrices, 

M  =  Mv  +  Mb  (21) 

C  =  Cv  +  Cb  (22) 

Combining  these  two  summations  with  the  total  damping  matrix,  D(u),  the  total  force/moment 
on  the  submarine  body  with  respect  to  the  body  frame  is, 

Mz>  +  C(v)v  +  D(u)u  +  g(rj)  +  g(v)ui  =  r  (23) 

Nevertheless,  the  analysis  of  the  motion  of  the  submarine  with  respect  to  the  earth-fixed 
frame  is  of  ultimate  interest,  so  use  the  translation  derived  in  equation  (4).  This  transfor¬ 


mation  provides  the  vehicular  dynamics  in  the  earth-fixed  frame, 

Mp(r))ij  +  Cv(rj,  v)i]  +  D^rj,  v)i]  +  gp(r))  +  g(u)ui  =  tv  (24) 

Mv(rj)  =  J-T(v)MtotJ-\rj)  (25) 

Cr,(vi  rj)  —  J~T(rj)  Ctotv  -  MtotJ-\V)j(r))  J_1(? 7)  (26) 

Dv(yV)  =  J-t(V)D(u)J-\V)  (27) 

9v(v)  =  J~T(v)d(v)  (28) 


where, 
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A. 2  Derivation  of  Control  Surface  Modeling 

Moments  are  produced  due  to  the  offset  of  the  actuation  rod  to  the  center  of  pressure  of  the 
lift  and  drag  forces.  By  assuming  a  constant  center  of  pressure,  the  moment  arm  for  the 
hydrodynamic  forces  also  remains  constant.  Finally,  the  coefficient  of  lift  for  small  angles  of 
attack  can  be  represented  by  the  Hoerner  empirical  fin  lift  coefficient  (ch),  [81], 

1  A, 

1.87T  7T  111 

where  As  is  the  surface  area  of  the  control  surface  and  hs  is  the  height  of  the  control  surface. 
Since  the  deflection  of  the  control  surface  is  restricted  to  ±30  degrees,  the  Hoerner  angle 
assumption  should  hold  in  this  case.  Altering  the  equation  of  lift  and  accounting  for  the 
duality  of  the  control  surfaces  leads  to, 

FL,m  —  pU2AsCs  (31) 

where  U  is  the  coupled  body-fixed  velocity  function  which  is  dependent  on  both  coupled 
velocities  as  well  as  control  surface  angles  of  deflection.  Expanding  the  coupled  body-fixed 
velocity  function  and  accounting  for  the  center  of  pressure  moment  arm,  the  hydrodynamic 
moment  equation  is  found, 


(30) 


Mf  —  pAsCslp  [u25s  +  uw  ±  lpuq\  (32) 

where  lp  is  the  moment  arm  to  the  center  of  pressure,  5S  is  the  control  surface  deflection,  u 
is  the  surge  velocity,  w  is  the  heave  velocity  and  q  is  the  pitch  velocity.  Coupled  with  the 
control  surface  dynamics,  the  hydrodynamic  moments  about  the  control  surface  produce  the 
required  actuation  torque  for  both  steady  state  and  transient  actions. 

An  interesting  note  on  the  linearized  hydrodynamic  moment  is  that  the  moment  function 
behaves  as  a  linear  compliance.  Combining  terms  to  form  the  linearized  moment  compliance 
parameter,  the  function  is  simplified  to  a  linear  impedance  element  as  shown  by  equation 
(33). 

Mfj  =  k09s  =>  kQ  =  pAgCglpU 2  (33) 


A. 2.1  Control  Surface  and  Electromechanical  System  Modeling 


This  control  surface  system  model  assumes  the  actuation  system  will  couple  into  rotational 
inertia,  Jm,  which  is  connected  to  the  control  surface  by  a  drive  shaft  having  stiffness, 
K.  On  each  end  of  the  drive  shaft,  bearing  losses  are  accounted  for  by  linear  frictional 
damping  coefficients  B i  and  B2.  The  control  surface  is  lumped  as  a  single  inertia  with 
moment  of  inertia,  Js.  Additionally,  the  linearized  ocean  compliance,  from  equation  (33)  is 
also  incorporated  into  the  two-port  model  as  another  applied  torque  on  the  control  surface 
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Control  Surface 
Angular  Velocity  -  cjcs 

-)  — 

Linearized  Hydrodynamic 
Compliance  -  K0 


Figure  22:  Simplified  layout  of  control  surface  system. 

R :  bi 
J 

Sf :  ub  | - ^  0  - - — ^  1  |- 

r  Jr. 

Se:Ta  l:Jms  Se:Td 


C  :  k/s  R :  b2  C  :  k0  Is 


Figure  23:  Bond  graph  form  of  the  control  surface  system  with  ideal  torque  actuator. 


inertia  [53]. 


A. 2. 2  Bond  Graph  Approach  for  Generic  Representation  of  Components 


Since  force  components  are  assumed  to  be  linear  including  the  ocean  hydrodynamic  moments, 
the  resulting  equation  is  a  linear  ordinary  differential  equation.  Figure  23  contains  the 
topology  of  the  energy  flow  in  the  control  surface  system  in  bond  graph  form.  The  state 
equations  for  the  control  surface  system  dynamics  are  derived  using  the  bond  graph  and  take 
the  form, 


hs 

1 

O 

k 

Os 

0 

0 

hm 

0\j  Jm 

—k 

1 

r-se 

1/Jm 

0 

- b2/Js 
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kQ 

1  /Js 
0 
0 


ha 

'  1 

0 

0  " 

'  Td  ' 

Os 

+ 

0 

0 

0 

hm 

0 

0 

1 

ojb 

T 

i 

0 

0 

0 

1  a 

(34) 


Due  to  the  linearity  of  the  control  surface  differential  equation  in  equation  (34),  the  differ¬ 
ential  equations  can  be  readily  expressed  as  a  transfer  function  relating  the  input  torque  to 
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Figure  24:  Magnitude  of  the  s-domain  representation  of  the  Pierson-Moskowitz  ocean  spec¬ 
trum. 


the  output  surface  deflection, 

Ldcs  __  S 

Ta  n4s4  +  n3s3  +  n2s2  +  n^s  —  ka 
where  the  denominator  coefficients  are  defined  by, 


n4 


n3 


n2 


n  i 


J$Jn 


k 

biJs  b'i  'fri 

~T  +  ~k~ 
bih 
k 


T  J s  Jm~f 


Jmko 


1  blk°  A  U 

bl~lT  +  b 


(35) 


(36) 


A. 3  Pierson-Moskowitz  Spectrum 


Many  different  models  for  ocean  wave  spectrum  exist,  but  for  the  purpose  of  this  derivation, 
the  Pierson-Moskowitz  ocean  spectrum  is  used  [65,  7].  Typically  the  Pierson-Moskowitz 
ocean  spectrum  is  given  by, 

=  0-0081^-3.109//^  (37) 

where  c o  is  the  frequency  in  radians  per  second,  g  is  the  acceleration  due  to  gravity,  and  hs  is 
the  significant  wave  height  in  meters.  However,  since  all  models  are  derived  in  the  s-domain, 
a  s-domain  representation  of  the  Pierson-Moskowitz  spectrum  is  required.  A  fourth  order 
s-domain  transfer  function  has  been  derived  in  [7]  and  is  expressed  as, 


G(s) 


_ KgS2 _ 

s4  +  cqs3  +  a2s2  +  a3s  +  a4 


(38) 


where  Kn  =  1.341,  a4  =  0.641,  a2  =  0.648,  a3  =  0.135,  and  a4  =  0.0641.  The  s-domain 
Pierson-Moskowitz  ocean  spectrum  magnitude  is  shown  in  Figure  24.  The  ocean  disturbance 
spectrum  contains  discernible  torques  from  0.18  rad/s  to  1.25  rad/s  with  a  peak  at  0.4  rad/s. 
Referring  back  to  Figure  11,  the  magnitude  of  the  natural  response  of  the  control  surface 
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Figure  25:  Model  geometries  used  by  Watt  [5]  to  identify  the  rotational  and  mass  moment 
derivatives. 


system  is  near  the  peak  in  this  band  of  frequencies. 


A. 4  Submarine  Parameter  List  for  Simulations 

The  parameters  used  in  this  report  are  based  on  the  geometries  of  the  submarine  model 
derived  by  Watt  [5].  However,  the  dimensions  used  throughout  this  report  are  based  on 
the  modern  Virginia  class  submarine  which  has  approximate  dimensions  given  by  Table  3. 
Parameter  scaling  allowed  for  the  calculation  of  new  mass  coefficients  and  moments  via  the 
existing  model  geometries  shown  in  Figure  25. 


Table  3:  Submarine  Model  Constant  Parameters  for  Virginia  Class  Vessel 


Parameter 

Value 

Description 

l 

110  m 

Length  of  submarine 

d 

1/ 8.75  m 

Hull  diameter  (semi-minor  axes) 

9 

9.81  m/s2 

Acceleration  of  gravity 

P 

1010  kg/m3 

Salt  water  fluid  density 

B  Synthesis  Techniques  Overview 


In  two-port  systems  effort,  e,  and  flow,  /,  variables  typically  represent  power  variable  pairs 
such  as  voltage-current,  force  velocity,  etc.  Consequently,  any  design  problems  including 
power  transmission  can  be  effectively  modeled  as  an  aggregate  of  two-port  subsystems  [43]. 
The  circuit  and  bond  graph  representations  of  a  two-port  are  illustratedi  in  Figure  26.  In 
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Table  4:  Hydrodynamic  and  added  mass  coefficients  for  submarine  dynamics  [5] 


Parameter 

Value 

Parameter 

Value 

Parameter 

Value 

Xg 

0.00 

y9 

0.00 

% 

0.55 

xb 

0.00 

Vb 

0.00 

Zb 

-0.55 

Ixx 

2.44  x  108 

lyy 

8.13  x  109 

Izz 

8.13  x  109 

Xu 

-3.69  x  105 

Y, 

-1.32  x  106 

Zu 

-1.91  x  103 

X. * 

-1.91  x  103 

Yp 

-2.16  x  107 

Zw 

-1.09  x  107 

*4 

6.98  x  104 

Y, 

1.52  x  106 

Zq 

-1.90  x  107 

Ky 

-2.16  x  107 

Mu 

6.98  x  105 

N, 

1.53  x  108 

Kp 

-2.18  x  108 

-1.90  x  107 

Np 

-1.67  x  108 

I\r 

-1.67  x  108 

Mq 

-7.49  x  109 

Nr 

-7.64  x  109 

Xuu 

1.00  x  104 

Y 

1  vv 

2.68  x  106 

ZyjW 

2.68  x  104 

y 

yvpp 

-2.08  x  106 

2.42  x  108 

Y 

yv  wp 

-1.65  x  105 

Y 

yvuq 

1.28  x  104 

Xqq 

2.45  x  108 

Xq\q\ 

3.70  x  104 

Y 

1  up 

3.36  x  106 

Yp\p\ 

-1.92  x  105 

Y 

1  ur 

7.43  x  106 

Yr\r\ 

2.29  x  108 

7 

^ wp 

-8.60  x  103 

7 

^ uq 

-8.84  x  106 

Zq\q\ 

-2.29  x  108 

Kup 

-3.13  x  107 

Xp\p\ 

-1.63  x  106 

XUr 

-1.02  x  107 

Kr\r\ 

-1.63  x  106 

Mpp 

2.11  x  107 

Mrr 

2.52  x  105 

Mwp 

-2.44  x  105 

Muq 

4.49  x  108 

Mq\q\ 

-1.06  x  1040 

Nup 

-4.73  x  107 

Xp\p\ 

-2.44  x  106 

N 

1  v  wp 

-1.83  x  106 

Nur 

-4.77  x  108 

Nr\r\ 

-1.06  x  1040 

XfisSs 

-1.27  x  105 

X$r$r 

-1.27  x  105 

Zss 

-1.58  x  105 

YSr 

1.46  x  105 

NSr 

-7.34  x  106 

Mss 

-7.34  x  106 

Xvr 

-n 

Xrp 

~Yt 

Xyjq 

Zw 

Y 

1  wr 

Xu, 

Y 

1  qr 

Xq 

Y 

1  wp 

Zw 

Y 

1  pq 

-Zq 

Zyjq 

-Xu, 

Zqq 

Xw 

Y 

y*-vp 

-Yu 

7 

ZJpp 

Yp 

7 

ZJrp 

Y, 

K 

±vwp 

~Yp 

K 

1^pq 

Kb 

K 

±vvq 

Yi-  +  Zq 

K 

J\yjr 

~Yr  ~  Zq 

K 

±\qr 

— Mq  +  Nr 

Mwq 

Xq 

Mvr 

Yp 

Mvp 

-Y, 

Mpr 

Kp  -  Nf. 

Nqr 

—  Kr 

Nvq 

1 

.  -Cr 

1 

Xpq 

—  Kp  +  Mq 

this  body  of  work,  the  impedance  synthesis  technique  relied  exclusively  on  two-port  systems. 


B.l  Canonical  Matrices 

Early  developers  of  n-port  matrix  theory  designated  certain  matrix  forms  as  “canonical”. 
A  detailed  derivation  of  two-port  matrix  forms  can  be  found  in  Huelsman  [43]  and  the 
matrix  forms  are  summarized  in  Figure  27.  The  synthesis  methodology  relies  heavily  on  the 
canonical  matrices  and  transformations  between  the  various  forms,  summarized  in  Figure 
28  [48].  Using  transformation  relationships,  the  system  matrix  can  be  defined  in  multiple 
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ei 

W 

Two-port 

fi 

Network 

6 

Effort  =  e  =  {voltage,  force,  torque,  pressure} 

Flow  =  f  =  {current,  velocity,  angular  velocity,  flowrate} 

(b) 


Figure  26:  (a)  Electrical  two-port,  (b)  Generalized  (multi-disciplinary)  two-port  with  bonds 
conveying  power,  P*  =  ej). 
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Figure  27:  Two-port  system  canonical  forms  with  descriptions  of  transfer  function  relation¬ 
ships. 
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Figure  28:  Relationships  between  canonical  two-port  systems.  Conversions  are  all  with 
respect  to  the  transmission  matrix,  which  can  be  easily  derived  from  the  system  bond  graph. 


ways  allowing  many  different  transfer  function  relationships  to  be  analyzed.  These  system 
techniques  have  been  extended  to  any  physical  domain  for  two-port  representation.  Due  to 
the  extension  to  multiple  physical  domains,  the  method  of  generating  the  two-port  system 
matrices  involves  the  use  of  bond  graph  topology. 


B.2  Impedance  Bond  Graphs 

Bond  graph  elements  represent  the  basic  linear  or  nonlinear  multi-disciplinary  elements  of 
realizable  systems.  For  those  unfamiliar  with  bond  graphs,  Karnopp  [8],  Redfield  [68,  69,  70] 
and  Connolly  [32,  17]  contain  overviews  of  the  subject.  Figure  29  presents  a  summary  of  the 
potential  bond  graph  elements  with  equivalent  impedances  based  on  the  causality  definitions 
of  the  bond.  For  example,  an  electrical  analog  model  for  the  ideal  torque  source,  Ta(t),  is 
shown  in  Figure  30(a).  Retrofitting  the  torque  source,  a  designer  can  remove  the  known 
ideal  torque  source  and  replace  it  with  an  unknown  impedance,  Za(s )  =  Ta/uia,  shown  in 
Figure  30(b)  as  electrical  analog  or  as  in  (c)  using  an  equivalent  bond  graph  form.  The 
synthesis  methodology  determines  the  unknown  impedance  function  in  terms  of  specified 
system  response  requirements. 
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Figure  29:  Basic  elements  of  bond  graph  notation  with  causal  definitions.  The  casual  defi¬ 
nition  implies  the  transfer  function  of  an  element. 
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Figure  30:  (a)  Electrical  analog  for  ideal  torque  actuation,  (b)  Electrical  analog  showing 
active  impedance  replacing  torque  source,  (c)  Bond  graph  form  of  system  in  (b). 
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Band-reject  frequency  response 


Figure  31:  Band  rejection  response  alteration  to  impose  on  the  natural  control  surface  system 
response. 


C  Derivation  of  Actuator  Synthesis 


The  detailed  synthesis  procedure  implemented  on  the  control  surface  model  is  discussed  in 
the  following  section.  This  process  draws  heavily  from  the  theory  discussed  in  Appendix  B. 


C.l  Specified  Frequency  Response  Dynamics  Derivation 


To  filter  out  ocean  disturbance  torques,  the  natural  response  of  the  system  is  modified  to 
have  a  band  reject  response.  This  modifying  function  is  shown  in  equation  (39). 


/ 

e 


A  +  1 
n 


s 

t^/ol 

s 

T2 


■  +  1' 

+T 


(39) 


where  T\  —  0.1  is  the  lower  limit  of  the  band  reject  function,  T2  =  5  is  the  upper  limit  of 
the  band  reject  function,  a  =  6  is  the  slope  factor,  n  —  2  is  the  filter  order.  This  function 
has  a  response  as  shown  in  Figure  31.  By  imposing  the  band  rejection  function  to  the 
natural  system  response,  the  desired  system  response  function  is  determined.  Combining 
these  functions  yields  equations  (40). 


dTF(s) 


1 

( Js  +  Jm)  s  —  +  B\  +  i?2 


—  + 1  A-  +  i 

OLT 1  1  Til  ot 

S-  +  1  '  S-  +  1 

n  t2 


(40) 


Merging  these  functions  produces  a  desired  transfer  function  that  will  not  radically  alter  the 
natural  response  of  the  system,  but  modifies  the  system  to  attenuate  disturbance  frequencies 
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Figure  32:  Comparison  of  the  desired  frequency  response  and  the  natural  system  response 
with  a  band  rejection  function. 


from  0.1  rad/s  to  5  rad/s.  The  desired  frequency  response  which  is  used  for  the  synthesis 
procedure  is  shown  in  Figure  32. 


C.2  Submarine  Actuator  Two-Port  Impedance 


Representation  of  the  control  surface  system  model  in  bond  graph  topology  is  shown  in 
Figure  33.  Using  the  bond  graph  representation,  the  two-port  matrices  of  each  junction 
are  formulated  using  transmission  matrix  form.  Beginning  with  the  right  side  of  the  bond 
graph  model,  the  matrices  are  constructed  for  each  element  and  parsed  into  series  matrix 


R :  bi  C :  k/s  R :  b2  C :  k0  Is 


Z(s)  l:JmS  Se:Td 


Figure  33:  Bond  graph  representation  of  control  surface  model  with  specified  unknown 
impedance  function. 
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Figure  34:  Comparison  of  the  desired  frequency  response  and  the  natural  system  response 
with  a  band  reject  function. 


multiplication, 
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(41) 


For  this  derivation,  the  drive  shaft  stiffness  K  is  allowed  to  approach  infinity.  This  assump¬ 
tion  simplifies  the  control  surface  model  by  eliminating  high  frequency  harmonics  from  the 
drive  shaft  compliance.  From  the  series  of  element  junctions,  the  impedance  function  that 
represents  the  entire  two-port  system  is  calculated.  The  total  function  is  represented  by 
equation  42. 


'  Ti  ' 

1  —  yr  +  (Ja  +  Jm)s  +  B\  +  B2 

'  Td  ' 

_  u  i  _ 

Ya(s)  Ya(s)[(Js  +  Jm)s-^  +  B1  +  B2]+l_ 

CS 

Now  that  both  the  two-port  model  of  the  system  and  the  desired  transfer  function  are 
specified  and  ensured  to  be  in  proper  form,  the  actuation  impedance  is  synthesized.  Due  to 
the  two-port  model  form,  we  have  four  equations  with  four  unknowns  to  solve.  The  elements 
H n,  ff12,  and  H2 1  are  equivalent  to  unspecified  unknown  frequency  responses,  and  H22  is 
equivalent  to  the  desired  frequency  response, 


Hn=>A 
H  \  2  =>  B 
H2\  =t  C 

H2  2  =>•  D 


Uu 
U12 
U2 1 


-+1 


(43) 


Solving  the  system  of  equations  (43)  provides  not  only  the  resulting  impedance  function  of  the 
unknown  impedance  function  Z(s),  but  also  the  other  relations  in  the  two-port  system  ( Uu , 
U 12;  U2 1).  The  resulting  frequency  response  for  the  transfer  functions  of  each  element  of  the 
two-port  matrix  are  summarized  in  Figure  34.  As  seen  in  the  graphs,  the  frequency  response 
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Frequency  Response  of  Actuator  Impedance  Function 


Frequency  (rad/s) 


Figure  35:  Frequency  response  of  the  synthesis  derived  actuation  impedance  function  for  the 
control  surface  system. 


for  H2 2  matches  the  specified  response.  Additionally,  the  resulting  impedance  function  for 
Z(s)  is  found  as  in  equation  (44).  This  impedance  function  relates  the  difference  in  velocity 
between  the  submarine  base  and  the  control  surface  to  a  torque  required  by  an  actuator 
system,  so  it  is  expressed  as  an  actuator  impedance, 


04s4  T  a3s3  T  a2s2  T  cqs  T  00 
M4  +  b3s3  +  b2s2  +  bis  +  60 


(44) 


where  a4  =  2.63  x  1018,  a3  =  3.42  x  1019,  a2  =  9.50  x  1019,  a,  =  4.55  x  1019,  a0  =  5.00  x  1018, 
b4  =  1.55  x  1013,  63  =  4.43  x  1013,  b2  =  4.72  x  1013,  61  =  2.22  x  1013,  and  b0  =  3.87  x  1012.  The 
actuation  impedance  function  has  a  frequency  response  as  seen  in  Figure  35.  The  actuation 
frequency  response  peaks  at  0.6  rad/s  which  means  that  the  torques  provided  by  the  actuation 
function  are  the  greatest  at  that  point.  A  positive  real  test  reveals  that  equation  (44)  is  a 
positive  real  function,  which  implies  a  purely  passive  realization  of  the  function  is  permissible. 
However,  further  investigation  indicates  that  positive  real  pole  removal  is  required  at  values 
other  than  s  =  0  or  s  =  00.  As  a  result  the  fourth  order  synthesized  impedance  requires 
either  complex  coupled  Brune  elements  or  active  components  for  realization  [28,  79]. 


C.3  Impedance  Expansion 

Using  partial  fraction  expansion,  the  actuation  impedance  function,  Za(s),  can  be  expanded 
into  separate  terms  so  that  individual  impedance  function  components  can  be  identified. 


Retrofit  Design  of  Submarine  Actuation  Systems 


Page  40 


I :  Ja  S  R:bi  C :  k/s  R:b2  C:k0/s 


Figure  36:  Synthesized  bond  graph  elements  separated  into  both  active  and  passive  compo¬ 
nents. 


Partial  fraction  expansion  of  equation  (44)  yields, 


•^a('S)  Z active^ T  Z passives) 

Z active  (^)  =  s2  s  f~  +  s2  s  ~  +  _s_  ,  _J_ 

k22  +  k21^~b21  fe32  +  fc31^  *>31  kil  641 

Zpassivei^)  ^11  4“  j  i  i~ 
fc12  *>12 

enabling  us  to  partition  the  actuator  impedance  into  separate  active  and  a  passive  compo¬ 
nents,  based  on  the  required  form  of  a  positive  real  function.  For  the  specific  case  examined 
here,  for  example,  we  can  obtain  passive  parameters  bn  =  1-70  x  105  N/s,  k\2  =  9.16  x  10' 
N-m/rad,  and  612  =  1.10  x  108  N/s,  and  active  parameters  k2\  =  4.82  x  106  N-m/rad, 
&21  =  1.61x10'  N/s,  k22  =  5.79x10®  N-m/rad,  k3 1  =  1.04x10'  N-m/rad,  631  =  2.49xl07lV/s, 
k32  =  1.73  x  107  N-m/rad,  k4 1  =  —8.98  x  107  N-m/rad,  and  641  =  —1.50  x  108  N/s. 


The  active  components  are  either  negative  impedances  or  non-realizable  complex  impedance 
elements  which  cannot  be  represented  by  simple  impedance  components.  The  passive  actu¬ 
ation  components  have  basic  impedance  relations  that  can  be  realized  as  physical  elements 
such  as  rotational  springs  and  dampers.  The  full  bond  graph  representation  of  the  derived 
impedance  function  within  the  framework  of  the  control  surface  system  is  shown  in  Figure 
36. 


C.4  Purely  Active  Realization  of  Synthesized  System 

If  the  passive  and  active  components  of  the  synthesized  impedance  function  are  not  separated, 
then  the  derived  impedance  function  is  ultimately  represented  by  a  single  modulated-effort 
source,  as  shown  in  Figure  37.  The  modulation  of  the  torque  source  is  prescribed  by  the 
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Figure  37:  Control  surface  system  realization  with  ideal  rotational  actuator. 


derived  actuation  impedance  function.  For  the  torque  source,  the  output  effort  (Ta)  is  related 
to  the  flow  signal  (wn)  via  the  transfer  function  of  the  active  elements.  Using  the  state 
equations  derived  for  the  active  elements  in  equation  46,  a  compensator  controller  function 
can  be  derived  in  state  space  form, 

0i2  =  can 
02i  =  can 
031  =  C^ll 
04i  =  can 

where  passive  parameters  are  bn  =  1.70  x  105  N/s,  fc12  =  9.16  x  10'  N-m/rad,  and  b 12  = 
1.10  x  108  N/s,  and  active  parameters  k2\  =  4.82  x  106  N-m/rad,  b2\  =  1.61  x  107  N/s, 
k22  =  5.79 x  106  N-m/rad,  k:i \  =  1.04 x  10'  N-m/rad,  i>31  =  2.49 x  10'  N/s,  k32  =  1-73 x  10'  N- 
m/rad,  hgi  =  —8.98  x  10'  N-m/rad,  and  641  =  —1.50  x  108  N/s.  Using  the  above  differential 
equations,  the  specified  torque  output  of  the  ideal  actuator  as  related  to  the  difference 
angular  velocity,  can  (difference  between  the  submarine  pitch  rate  and  the  control  surface 
rate),  results  from  the  summation  of  the  impedances  given  by, 

Ta  =  biiOJu  +  ki20i2  +  k2i02i  +  k-31031  +  kii04i  (47) 

Representing  equation  (47)  in  block  diagram  form  results  in  a  controller  algorithm  that 
utilizes  the  parameters  derived  directly  from  the  synthesis  method.  Figure  38  depicts  the 
controller  form  integrated  with  the  bond  graph  structure. 


ki2  n 

bRUl2 
*210. 
021 
kzia 
b3 1  : 
fc-41  a 

t>41 


21 


31 


fc2i  a 

_  fc3i  a 

k32  ' 


41 


31 


(46) 


C.5  Separated  Active  and  Passive  Subsystem  Realization 

As  in  the  active  suspension  example,  separation  of  the  passive  elements  reduces  the  control 
algorithm  by  eliminating  passive  elements  from  the  function  and  physically  realizing  them 
as  passive  components.  For  the  sake  of  comparison,  the  energy  consumption  should  be 
compared  between  both  the  pure  active  case  and  the  separated  passive  case. 
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Figure  38:  Torque  source  controller  algorithm  as  determined  by  the  synthesis  expansion 
incorporated  in  the  bond  graph  structure. 
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Figure  39:  Control  surface  system  realization  with  passive  components  separated. 


Flow-induced 
Torque  Disturbances 


Figure  40:  Block  diagram  representation  of  control  surface  system  and  compensator. 


By  separating  the  active  and  passive  components,  a  modulated  effort  source  and  a  passive 
spring  damper  system  represent  the  derived  impedance  functions,  as  shown  in  schematically 
in  Figure  39,  or  in  block  diagram  form  in  Figure  40.  The  modulation  of  the  torque  source  is 
governed  by  the  derived  actuation  impedance  function. 

For  the  modified  torque  source,  the  transfer  function  of  the  active  elements  without  the 
passive  components  relates  the  output  effort  (Ta)  via  the  flow  signal  (cun).  Using  the  state 
equations  derived  from  the  active  elements  in  equation  (48),  the  compensator  controller 
function  state  equations  are, 


y = wn  -  -  fc92i 

031  =  ‘•'ll  -  gr»31  -  fc«31  (48) 

041  =  kfLi  —  fjj-041 

where  the  parameters  remain  identical  to  those  given  above  for  equation  (46).  Again,  the 
specified  torque  output  of  the  ideal  actuator  is  related  to  the  difference  angular  velocity,  Wn, 
and  is  the  summation  of  the  impedances, 


Ta  —  bnUn  +  &21021  +  &31031  +  ^41041 


(49) 
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Representing  the  modified  controller  function  equation  (49)  in  block  diagram  form  results 
in  a  modified  controller  algorithm,  this  utilizes  the  parameters  derived  from  the  synthesis 
method.  Figure  41  depicts  the  controller  integrated  within  the  bond  graph  structure. 


C.6  Model  Reduction  of  Synthesized  Impedance 


The  fourth  order  complexity  of  the  synthesized  impedance  function  complicates  the  realiza¬ 
tion  procedure.  Although  the  complex  impedance  function  yields  positive  real  characteristics, 
a  purely  passive  realization  of  the  function  could  not  be  completely  derived  due  to  function 
complexity.  A  purely  passive  realization  of  the  fourth  order  system  with  complex  roots 
requires  coupled  passive  elements  as  shown  in  Temes  [79].  Coupled  passive  elements  are 
physically  difficult  to  realize  in  many  physical  domains  and  should  be  avoided  while  using 
synthesis  techniques. 


The  order  of  a  given  system  can  always  be  reduced  via  balanced  realization  with  increasing 
error  depending  on  the  magnitude  of  the  state  contributions.  The  interested  reader  can 
find  additional  coverage  of  the  balanced  realization  procedure  in  Appendix  D.  The  original 
actuator  impedance  function, 


.  a4s4  +  a3s3  +  a,2S 2  +  cqs  +  Oo 
b4s4  +  b3s3  +  b2s2  +  bis  +  b0 

can  be  converted  to  the  state  space  form, 


'  -2.867 

-3.055 

-1.433 

-0.25 

1 

1 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

1 

0 

0 

1.725e6 

5.628e6 

2.698e6 

2.873e5 

1.702e5 

(50) 


(51) 


Computing  the  balanced  realization  of  the  impedance  space  state  function  produces  both 
the  balanced  state  space  realization  as  well  as  the  error  G  vector.  The  error  G  vector  for 
Zjs)  is  given  by, 

G error  =  [  1-5235  0.8543  0.0972  0.0026  ]T 

The  order  of  magnitude  drop  between  the  second  and  third  entries  in  the  error  G  vector 
implies  that  an  accurate  realization  of  a  second  order  system  exists  for  Za(s).  The  relative 
contributions  from  the  balance  third  and  fourth  states  are  an  order  of  magnitude  less  than 
the  contributions  of  the  first  and  second  states.  The  reduced  order  Za(s)  system  is  procured 
by  elimination  of  the  low  contribution  states  shown  in  equation  (52). 


"  -0.7043 

-0.6297 

-1465 

Zred\S) 

0.6297 

-0.0484 

287.7 

-1465 

-287.7 

1.702e5 

(52) 


Conversion  of  the  second  order  state  space  representation  to  a  SISO  transfer  function  pro- 
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Figure  41:  Torque  source  controller  algorithm  with  passive  components  separated. 
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Figure  42:  Reduced  order  impedance  function  in  bond  graph  form  yielding  a  purely  passive 
system. 


duces  the  reduced  impedance  function, 

1.702  x  105s2  +  2.192  x  106s  +  6.498  x  105 
s2  +  0.7527s  +  0.4306 

The  reduced  actuation  function  has  a  frequency  response  as  shown  in  Figure  18.  For  compar¬ 
ison,  the  frequency  response  plot  also  contains  the  original  synthesized  impedance  function. 
Since  the  original  impedance  function  is  a  positive  real  function,  the  reduced  order  system 
also  is  positive  real  such  that  a  purely  passive  realization  of  the  function  must  exist.  Due  to 
the  order  reduction,  the  realization  of  the  impedance  function  with  desired  response  yields 
a  less  complicated  result  which  can  more  readily  be  physically  constructed.  Continuous 
fraction  expansion  of  the  Zred(s )  function  produces  equation  (54). 

Zred{s )  =  bn  -I — — — I — — - 1 -  (54) 

k\2  b±2  772215+621 


where  bn  =  1-70  x  105,  kn  =  2.06  x  106,  bn  =  4.36  x  106,  m2 1  =  6.917  x  106,  and 
621  =  1.932  x  106.  Since  all  element  coefficients  from  the  resulting  function  expansion  are 
positive  real  values,  the  total  function  is  purely  passive.  The  full  bond  graph  representation 
of  the  reduced  impedance  function  within  the  framework  of  the  control  surface  system  is 
shown  in  Figure  42. 
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Figure  43:  Original  pure  mechanical  domain  system  and  the  new  partially  electrical  domain 
system  with  unknown  equivalent  two-port  impedance  function. 


C.7  Equivalent  Impedance  Transformation 


To  find  the  equivalent  impedance  of  the  new  modified  bond  graph  via  secondary  synthesis,  the 
original  bond  graph  transmission  matrix  must  equate  to  the  modified  transmission  matrix. 
Figure  43  shows  the  original  bond  graph  with  the  flow  source  reverted  to  the  original  negative 
impedance  function  that  was  derived  and  the  pure  hydraulic  form  bond  graph.  Converting 
the  compensator  algorithm  back  into  an  impedance  function  ensures  that  any  necessary 
modifications  to  the  controller  are  accounted  for. 


For  each  bond  graph,  the  two-port  transmission  matrices  must  be  found  for  modified  portions 
of  the  bond  graph.  Since  no  changes  occur  beyond  the  synthesized  impedance  function, 
the  remaining  portions  of  the  bond  graph  are  not  needed  for  the  two-port  representation. 
Equation  (55)  illustrates  the  matrix  function  for  the  original  case  with  a  placeholder  gyrator 
element  added,  and  equation  (56)  displays  the  matrix  function  for  the  electrical  conversion. 


r  0 

Kt  1 

r  1 

0  1 
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-  ^red(^) 
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(55) 


In  the  case  for  the  electrical  conversion,  the  equivalent  two-port  impedance  is  unknown  and 
thusly  represented  as  a  matrix  of  unknowns. 


Z\2 

Z22 


0 

1 

Kt 


Kt 
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IS  KtZ  ill 
IS  K Tzn. 


(56) 


For  each  two-port  representation  to  have  equivalent  impedance  matrices,  the  derived  two- 
port  matrices  must  be  identically  equal.  Solving  for  the  unknown  Z  matrix,  the  equivalent 
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structure  is  found  in  equation  (57). 


r 

1 

I< f. 

^12 

— 

^red(^) 

Z2i 

z22 

0 

1 

(57) 


The  impedance  function  derived  in  the  above  equation  has  the  structure  of  a  1-junction 
where  the  impedance  is  composed  of  the  inverse  of  the  elements  off  the  junction.  However, 
due  to  the  gyrator  modulus,  there  is  an  additional  scaling  factor  on  all  of  the  impedance 
elements.  To  compensate  for  the  gyrator  modulus,  the  coefficients  of  the  impedance  elements 
must  be  scaled  accordingly.  The  resulting  equivalent  impedance  function  takes  the  form, 

Ze,(s)  =  ^ -  (58) 

Kt  +  m  12  K'tS  +  bl2K2T  +  s/k2iKl\s/b2iKl 

where  bu  =  1-70  x  105,  m12  =  2.06  x  106,  bV2  =  4.36  x  106,  k21  =  6.917  x  106,  b21  =  1.932  x  106 
and  Kt  is  any  positive  real  number.  With  these  modified  parameters  and  the  structure  from 
the  two-port  secondary  synthesis,  the  new  bond  graph  structure  can  be  generated. 

Figure  44  shows  the  new  bond  graph  layout  of  the  electrical  conversion  realization  of  the 
control  surface  ocean  disturbance  rejection  compensator.  Realization  of  the  gyrator  and 
electrical  domain  portion  of  the  bond  graph  yields  Figure  21.  The  gyrator  represents  the  ideal 
torque  source  actuation.  Moving  the  compensation  function  from  the  mechanical  domain  to 
the  electrical  domain  produces  an  internal  voltage  filter  which  indirectly  filters  the  torque 
supplied  to  the  control  surface  system. 


D  Balanced  Realization 


In  the  case  of  several  of  the  derived  impedance  functions  found  throughout  this  report,  the 
order  of  the  resulting  function  was  the  magnitude  of  both  the  order  of  the  model  function 
and  the  desired  response  function.  For  example,  a  second  order  model  with  an  unknown 
impedance  function  that  is  run  through  the  synthesis  procedure  with  a  desired  response  of 
a  second  order  system  procures  a  fourth  order  impedance  function.  Higher  order  impedance 
functions  tend  to  be  physically  difficult  to  realize  and  even  more  difficult  to  physically  man¬ 
ufacture.  Additionally,  high  order  functions  that  are  found  to  be  positive  real  often  cannot 
easily  be  decomposed  into  purely  passive  elements  without  significant  function  manipulation. 

One  potential  method  for  impedance  model  simplification  is  order  reduction  of  the  function 
[58].  Reliable  reduction  of  model  order  can  be  achieved  by  internally  balancing  the  model 
system.  This  balanced  realization  contains  the  original  number  of  states,  but  readjusts  the 
state  space  values  so  that  the  relative  contributions  from  each  state  are  in  descending  order. 
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Figure  44:  Passive  compensation  as  realized  in  electrical  form  inside  electromechanical  ac¬ 
tuator. 
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Consider  a  large  order  impedance  transfer  function  Z(s)  which  must  be  either  proper  or 
strictly  proper.  By  the  definition  of  the  transfer  function,  the  function  Z (s)  can  be  converted 
from  transfer  function  form  to  the  state  space  form. 

ry,  .  _  N(s)  _ ^  x  —  Ax  +  Bu 
[S)~~D{s)^  y  =  Cx  +  Du 

A  state  coordinate  transformation  (x  =  Tx)  can  be  performed  on  the  state  space  represen¬ 
tation  of  the  function  Z(s).  This  similarity  transformation  results  in  an  equivalent  system 
as  shown. 

'x  =  TAT-'x  +  TBu 
y  =  CT~lx  +  Du 

In  the  case  of  a  balanced  realization,  the  unique  transformation  T,  when  applied  to  both  the 
observability  gramians  (L0)  and  controllability  gramians  (Lc),  results  in  modified  gramians 
that  are  equivalent.  The  modified  gramians  are  given  as, 

L0  =  T~tL0T -1 
Lc  =  TLcTr 

In  the  case  of  an  internally  balanced  realization  when  the  modified  gramians  are  equivalent, 
the  resulting  structure  of  the  modified  gramians  is  in  diagonalized  form, 

L o  Lc  diuy^c q,  cj2i  ■  •  •  9n) 

where  the  diagonal  terms  are  in  descending  order  (g\  >  g2  >?  >  gn).  Reduction  of  system 
order  eliminates  states  from  the  balanced  realization  system  according  to  the  diagonalized 
gramians  form.  A  proper  model  reduction  identifies  an  order  of  magnitude  drop  in  the 
diagonalized  matrix  as  an  indication  of  decreasing  state  contributions.  The  states  of  the 
balanced  realization  at  this  division  are  eliminated,  which  reduces  the  model  order.  A 
reduction  to  any  order  is  possible,  however  reduction  errors  become  ever  more  present  with 
increasing  order  reduction, 


||G(s)  —  Gapproa;(s)||oo  —  2  ^  Afc 

k=m-\- 1 

The  errors  of  the  model  reduction  are  bounded  by  the  sum  of  the  eliminated  states  from  the 
diagonalized  matrix  gramian  forms. 


D.2  Method  to  Determine  Balanced  Realization 

To  determine  the  unique  transformation  T  that  yields  the  balanced  realization,  an  algorith¬ 
mic  process  that  produces  a  decomposition  of  one  of  the  gramians  is  found  to  produce  equiv¬ 
alent  controllability  and  observability  gramians  [58].  The  following  procedure  for  procuring 
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the  unique  transformation  is  implementable  in  the  MATLAB  environment  with  the  “balreal” 
command. 


1.  The  observability  gramian  must  be  factored  into  two  equivalent  matrices  with  one 
being  orthogonal  to  the  other. 

L0  =  RtR 

where  the  matrix  R  cannot  be  singular. 

2.  Combining  the  R  matrix  found  from  the  observability  gramian  with  the  controllability 
gramian  produces  an  altered  controllability  matrix.  Decomposing  the  new  matrix  into 
Eigenvalue  and  eigenvector  matrices  yields  the  following: 

RLcRt  =  UA2Ut 

where  U  are  the  eigenvectors  which  combine  to  form  an  orthogonal  unitary  matrix  and 
A  is  a  diagonal  matrix  comprised  of  the  Eigenvalues. 

3.  The  following  procedure  yields  the  unique  transformation  that  produces  the  balanced 
realization.  The  matrix  T  is  defined  as  follows: 

T  =  R~lU  A1' /2 


E  Passivity  Conditions 


Conditions  exist  which  sufficiently  indicate  that  a  given  impedance  function  is  realizable 
via  strictly  passive  elements  [28].  However,  note  that  satisfying  these  conditions  does  not 
imply  that  the  passive  realization  of  the  system  is  succinctly  derived  with  straightforward 
techniques.  Complex  manipulations  and  couplings  of  impedance  functions  are  sometimes  re¬ 
quired  to  generate  a  purely  passive  realization,  especially  in  the  case  of  high  order  impedance 
functions.  The  following  definitions  provide  the  outlines  of  the  composition  of  passively  re¬ 
alizable  systems  in  terms  of  all  bond  graph  elements. 


E.l  One-Port  Passive  System  Conditions 

The  location  of  the  poles  and/or  zeros  in  a  given  polynomial  function  is  an  important  concept 
when  dealing  with  the  complex  s-domain.  Many  insights  into  the  function  response  and  type 
derived  from  the  location  of  poles  and  zeros  in  the  complex  s-domain. 


Definition  1  A  given  polynomial  function,  A(s),  in  the  s-domain  is  said  to  be  a  Hurwitz 
polynomial  if 

A(s)  —  oq  ~t~  ois  T  . . .  T  ans 
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1.  A(s)  is  a  real  polynomial,  meaning  that  all  values  at  are  real  coefficients. 

2.  All  of  the  zeros  of  A(s)  are  in  the  left  half-plane.  Any  zeros  that  lie  on  the  imaginary 
axis  must  be  simple  zeros.  (Re(s)  <  0). 


The  polynomial  A(s)  is  said  to  be  a  strictly  Hurwitz  polynomial  if  all  the  zeros  are  in  the  left 
half  plane  (LHP)  without  any  zeros  which  lie  on  the  j- axis.  (Re(s)  <  0). 


Definition  2  A  one-port  impedance  function,  Z(s),  is  said  to  be  passively  realizable  if  the 
given  impedance  function  is  found  to  be  a  positive  real  function  [28]. 


1.  The  impedance  function,  Z(s),  must  be  a  real  rational  function  of  s. 

y(  \  —  a°  +  aiS  +  •  •  •  +  anSn 

bo  +  b±s  +  . . .  +  bmsm 

where  the  coefficients  ai  and  bj  are  all  real  values. 

2.  For  all  real  values  of  u,  the  following  condition  is  satisfied: 

R eZ(juj)  >  0 

3.  All  the  poles  of  Z(s)  are  in  the  closed  LHP  of  the  s-plane. 


The  third  condition  of  positive  realness  on  the  pole  locations  can  also  be  expressed  in  terms 
of  Hurwitz  polynomials.  Depending  on  the  polynomial  type,  the  third  condition  is  satisfied  if 
the  following  conditions  are  met: 


1.  The  denominator  of  Z(s)  must  be  either  a  Hurwitz  or  strictly  Hurwitz  polynomial. 

2.  If  the  denominator  polynomial  is  only  Hurwitz,  then  the  zeros  on  the  jiv-axis  must  be 
simple.  Additio7ially ,  the  poles  as  s  — »  oo  must  also  be  simple. 


E.2  Methods  for  Determination  of  Positive-Realness 

Simple  tests  have  been  found  which  provide  insight  into  the  positive  real  nature  of  a  given 
function  [79].  The  first  condition  of  positive  realness  is  easily  determined  via  inspection  and 
verification  of  each  polynomial  coefficient  being  real  valued.  The  second  and  third  conditions 
are  more  difficult  to  ascertain  and  cannot  be  immediately  found  via  inspection. 

A  technique  which  is  a  sufficient  condition  for  proving  Re(jcj)  >  0  (condition  2  for  positive 
real  functions  given  above),  is  Strum’s  theorem.  Strum’s  theorem  provides  a  method  through 
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which  the  number  of  sign  changes,  meaning  a  simple  real  zero  value,  occur  on  the  positive 
axis.  Strum’s  theorem  presents  the  following  condition: 

P(x)  =  anxn  +  an- ixn_1  +  •  •  •  +  a±x  +  ao  >  0 

where  x  >  0  and  the  function  P(x)  is  found  by  division  of  the  function  of  interest  into  even 
and  odd  polynomial  values  for  both  the  numerator  and  denominator.  The  terms  are  then 
combined  together  as  follows  to  produce  P(x)  where  m2  =  x: 

P(uj2)  =  Neven(ju)D  even  {ju)  -  N odd(j  in)  Dodd(j  in) 

The  condition  of  Strum’s  theorem  requires  that  P(x)  contain  no  simple  real  zeros  on  the 
positive  interval,  which  satisfies  Re  (jin)  >  0. 


Sufficient  conditions  for  the  satisfaction  of  the  third  condition  of  positive  realness  can  also 
be  found  via  the  manipulation  of  numerator  and  denominator  of  the  function  of  interest. 
The  following  modification  of  the  function  Z(s)  must  have  all  zeros  in  the  open  LHP: 


mod 


(s) 


N(s)  N(s)  +  D(s) 

D(s)  +  D(s) 


The  modified  function  above  implies  that  N(s)  +  D(s)  must  be  strictly  Hurwitz  in  order  to 
have  all  zeros  in  the  open  LHP.  Therefore,  the  addition  of  the  numerator  and  denominator 
provides  a  quick  test  of  the  third  condition  of  positive  realness  of  a  given  function. 


E.3  Extension  to  Two-port  Systems 

The  prior  discussion  dealt  only  with  the  possibility  of  passive  realization  for  a  given  one- 
port  impedance  transfer  function.  The  following  definition  extends  the  previous  definitions 
to  encapsulate  two-port  functions  [79].  As  in  the  case  of  the  one-port  system,  a  two-port 
system  must  be  a  positive  real  function  for  a  purely  passive  realization  to  exist. 


Definition  3  A  two-port  impedance  function  system  given  by  a  two-port  matrix  Z(s)  is 
positive  real  if  the  following  conditions  are  satisfied: 


1.  Z(s)  exists  and  is  analytic  on  the  interval  of  a  >  0,  where  a  is  the  real  portion  of  s 
( s  =  o  +  ju). 

2.  Z*(s)  =  Z(s*)  on  the  interval  a  >  0,  where  the  superscript  asterisk  indicates  the 
complex  conjugate. 

3.  ZH(s)  >  0  on  the  interval  a  >  0,  where  ZH(s)  designates  the  Hermitian  portion  of  the 
matrix. 


The  above  conditions  for  positive  realness  of  two-port  systems  can  also  be  expressed  in  a 
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perhaps  more  readily  understood  form.  The  following  conditions  are  equivalent  to  the  initial 
three  of  the  definition: 

1.  All  elements  of  Z(s)  are  real  rational  functions. 

2.  Z(s)  contains  no  poles  in  the  range  of  a  >  0. 

3.  The  poles  of  Z(s)  on  the  jco-axis  are  simple  values. 

4-  For  the  poles  on  the  joj-axis,  the  residues  of  each  element  must  be  simple  values. 

5.  ZH(juj)  >  0  for  all  real  values  of  u. 
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